202103.00039v1 


chinaXiv 


干旱 


ChinaXiv 合 作 期 刊 


研 究 


ARID ZONE RESEARCH 


基于 InVEST 模 型 的 1999 一 2016 年 麻 塔 流域 
碳 储量 变化 及 空间 格局 研究 


法 12 
JE ’ 


Kw, BO H R” 


(1. 西北 农林 科技 大 学 林学 院 , 陕 西 杨凌 712100; 2. 西北 农林 科技 大 学 黄土 高 原 土 壤 侵 蚀 与 旱地 农业 


国家 重点 实验 室 ,陕西 杨凌 712100; 3. 中 国 科学 院 水 利 部 水 土 保持 研究 所 ,陕西 杨凌 712100) 


摘 要 : 区 域 碳 储量 是 生态 系统 功能 的 重要 度量 指标 ,探索 土地 改造 对 区 域 碳 储量 变化 的 影响 ,对 于 协调 区 域 生 
态 建设 和 生态 产业 发 展 具有 重要 意义 。 基 于 InVEST 模 型 和 地 理 信息 系统 技术 ,研究 了 陕西 省 延安 南部 麻 塔 流域 


1999 一 2016 年 土地 结构 改造 过 程 中 区 域 碳 储 量变 化 ,并 探讨 坡度 、 坡 向 、 坡 位 对 碳 储量 空间 分 布 的 影响 。 结 果 表 
明 : 麻 塔 流域 18 a 的 土地 改造 使 得 区 域 碳 储 量 增加 1688.36 Mg( 碳 密度 增加 6.92 Mg hm?) ,总 固 碳 功能 提升 约 


7.63%。 和 森林 草地 和 果园 土地 类 型 面积 的 增加 是 植被 转换 后 流域 景观 碳 储量 提升 的 主要 贡献 者 。 碳 储量 增加 的 
空间 位 置 主要 分 布 在 半 阴 坡 .中 上 坡 位 以 及 坡度 10"~30?。 本 研究 认为 林 草 植被 建设 和 经 济 果 林 建 设 两 者 都 有 利 
于 增强 麻 塔 流域 景观 固 碳 能 力 。 麻 塔 流域 土地 转换 模式 能 够 协调 区 域 生态 建设 和 生态 产业 的 发 展 ,在 黄土 丘陵 区 


具有 较 好 的 推广 价值 。 
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生态 系统 固 碳 作为 最 重要 的 生态 系统 服务 功 
能 之 一 ,对 于 降低 大 气 中 CO0, 等 温室 气体 浓度 .调节 
区 域 小 气候 ,减缓 全 球 气 候 变 化 以 及 维持 生态 平衡 
具有 至 关 重 要 的 作用 。 近 年 来 ,对 于 生态 系统 碳 
储量 及 空间 格局 的 研究 ,成 为 各 国政 府 与 国内 外 专 
家 关注 的 焦点 。 因 此 ,研究 区 域 矶 储量 及 其 时 空格 
局 的 变化 对 于 自然 资源 管理 和 生态 决策 具有 重要 
的 参考 价值 与 指导 意义 。 

目前 , 碳 储量 估算 方法 多 种 多 样 , 且 具 有 尺度 
差异 性 ”。 对 于 中 小 尺度 碳 储量 估算 大 多 采用 野外 
调查 仪器 测量 .数理 统计 等 方法 5 。 但 这 类 方法 
仅 适 用 于 中 小 尺度 碳 储量 估算 ,费时 费力 , 且 研 究 
结果 呈 静 态 ,无 法 准确 反映 研究 区 碳 储量 的 动态 变 
化 及 空间 格局 "”。 大 尺度 的 碳 储量 估算 ,一 般 基 于 
RS 和 GIS 技 术 建 立 模 型 来 研究 InVEST 模 型 是 
一 种 可 量化 多 种 生态 系统 服务 功能 的 综合 评 佑 模 
型 ,已 广泛 用 于 多 个 地 区 ,如 美国 ,亚洲 印度 
尼 西 亚 " 和 非洲 坦桑尼亚 等 ,及 中 国 汶川 地 震 
灾区 "中 长 江 中 下 游 平 原 " 下 北京 山区 "祁连山 


收 稿 日 期 : 2020-04-07; 修订 日 期 : 2020-08-05 


保护 区 :等 区 域 。 因 此 ,将 AreGIS 与 mVEST 模 型 
相 结合 ,可 更 快速 .简便 .准确 地 估算 区 域 碳 储量 及 
空间 动态 变化 ,并 以 地 图 形式 直观 展现 出 来 。 

黄土 丘陵 区 是 我 国 水 土 流失 严重 的 区 域 ,也 是 
生态 环境 脆弱 的 区 域 之 一 。20 世 纪 70 年 代 以 来 的 
水 土 保 持 工程 与 1999 年 以 来 的 退耕 还 林 工 程 , 极 大 
地 提升 了 该 地 区 的 植被 覆盖 度 。 国 内 外 学 者 对 黄 
士 丘陵 区 植被 固 碳 功 能 及 其 变化 进行 了 大 量 研究 ， 
这 些 研究 多 数 是 聚焦 于 水 土 保 持 工程 以 及 退耕 还 
林 工 程 的 固 碳 效益 。 然 而 ,在 当前 新 农村 建设 与 
生态 农业 等 政策 的 驱动 下 ,黄土 丘陵 区 土地 利用 政 
策 发 生 了 明显 的 变化 ,过 去 种 植 的 水 土 保 持 林 正在 
BOER (Malus pumila) 等 经 济 林 替 代 , 这 种 现象 在 延 
安 周 边 地 区 尤其 突出 到 。 因 此 ,探索 如 何 量化 经 济 
果林 建设 对 区 域 回 碳 能 力 的 影响 是 当前 亚 待 解决 
的 一 个 科学 问题 。 

本 文 以 黄土 丘陵 区 典型 代表 一 一 延安 麻 塔 流 
域 为 例 , 该 流域 从 1999 一 2016 年 ,苹果 种 植 面积 持 
续 增 加 ,在 当前 土地 利用 政策 驱动 下 ,为 植被 恢复 
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方式 的 转变 对 区 域 固 碳 功能 影响 的 研究 提供 了 一 
个 理想 的 对 象 。 因 此 ,本 文 试 图 通过 InVEST 模 型 研 
究 该 流域 1999 一 2016 年 生态 系统 固 碳 能 力 及 其 在 
空间 上 的 分 配 特征 ,探索 流域 生态 系统 固 碳 能 力 变 
化 及 其 分 布 格局 ,以 期 为 黄土 丘陵 区 植被 建设 和 生 
态 产 业 的 可 持续 发 展 提供 科学 依据 。 


1 研究 区 概况 与 研究 方法 


1.1 研究 区 概况 

麻 塔 流域 (109.5268°N ,36.4695°E) 位 于 陕西 省 
延安 市 宝塔 区 南部 (图 1) ,为 典型 的 黄土 丘陵 地 
带 。 流 域 面积 2.4 km ,海拔 1185~1379 m, BER 
热 ,冬季 寒冷 干燥 ,为 干旱 半 干 旱 气 修 , 洪 在 蒸 散 为 
544 mm, 年 均 气 温 在 8.8 左右 ,年 均 降 雨量 523 
mm ,全 年 降水 主要 集中 在 6 一 9 月。 该 流域 土壤 


mm 林地 ml 果园 访 计 
灌木 本 农田 d 
EE 草地 国 裸 地 


图 1 麻 塔 流域 位 置 及 1999 年 和 2016 年 土地 利用 类 型 
Fig. 1 Location of Mata Watershed and its land use types 


类 型 为 黄 绵 土 ,主导 产业 为 旱地 苹果 ,为 当地 居民 
支柱 产业 。 
1.2 土地 利用 数据 来 源 

麻 塔 流域 土地 利用 图 (1999 年 和 2016 年 ) 是 通 
过 解 译 美 国 陆地 卫星 1999 年 8 月 Landsat 7 和 2016 
年 9 月 Landsat 8 影像 的 遥感 数据 获得 。 遥 感 影像 选 
择 Landsat 7 影像 的 1~6 波 段 ;选择 Landsat 8 影像 的 
2~7 波 段 。 获 取 的 遥感 波段 经 过 地 形 校正 、. 辆 射 定 
标 、 辐 射 校正 ,将 它们 的 灰 度 值 转换 为 地 表 反 射 率 
数值 。2016 年 的 土地 利用 图 是 通过 野外 调查 195 个 
地 面 点 数据 和 处理 后 Landsat 8 影像 的 2~7 波 段 进行 
监督 分 类 (随机 森林 算法 ) 完 成 的 ,最 终 的 预测 准确 
率 为 74%(Kappa 值 为 0.68)。1999 年 十 地 利用 图 是 
将 上 述 训 练 的 随机 森林 分 类 模型 应 用 到 1999 年 
Landsat 7 影像 1~6 波 段 上 得 到 的 预测 结果 。 土 地 类 
型 包括 6 类 :林地 灌木 .草地 .果园 .农田 、 其 他 ( 道 
路 + 房屋 + 裸 地 )。 

由 表 1 结 果 显 示 , 森 林 主 要 转化 为 灌木 草地 和 
果园 ;灌木 主要 转化 为 草地 和 果园 ;草地 农田 及 其 
他 主要 转化 为 果园 。 在 6 种 土地 利用 类 型 中 ,森林 、 
草地 .果园 的 转 和 面积 大 于 转 出 面积 。 而 灌木 K 
FA .其 他 的 转 入 面积 小 于 转 出 面积 。1999 年 和 2016 
年 ,整个 流域 6 种 土地 利用 类 型 中 均 为 果园 面积 最 
大 ,1999 年 占 总 体面 积 的 34.7% ,2016 年 占 总 体面 
积 的 47.4%。 

13 流域 碳 储量 及 其 估算 

陆地 生态 系统 碳 储量 主要 包括 地 上 生物 碳 储 
E (Caoc) x HE F Æ W UR i E (Ce) .土壤 碳 储量 
(Coi) JECA HURREE (Caua) ORBA HEP Catone HH 
H EA FP Yi HEL A PET D SE H HS k 
HE 3 Cuw 指 植物 活性 根系 中 的 碳 储量 ; Cw 指 矿质 土 


表 1 1999—2016 年 麻 塔 流域 土地 利用 类 型 面积 转移 矩阵 


Tab.1 Area transfer matrix of land use types in Mata Watershed from 1999 to 2016 /hm2 
1999 年 cau = 
森林 灌木 草地 果园 KH 其 他 合计 
森林 5.21 2.07 1.96 1.94 0.27 0.15 11.60 
灌木 7.53 15.23 17.96 26.32 3.08 6.66 76.78 
草地 0.83 1.06 4.99 10.37 1.55 2.66 21.46 
果园 6.10 8.59 20.43 43.07 2.68 3.27 84.14 
农田 0.11 0.53 3.38 12.01 1.28 1.85 19.16 
其 他 0.15 0.35 5.86 21.26 1.05 0.67 29.34 
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114.97 9.91 15.26 242.48 
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壤 和 有 机 土壤 中 的 碳 储 量 ; Cu 包括 死亡 植被 或 凋 
落 物 中 的 碳 储量 。 碳 储量 计算 公式 为 : 

Cinai = Case FC peia F Gaai FC aag 

SUP : Coa RIR BBA E (Mg) ;Cw 表示 地 上 生物 碳 
fifi (Mg) ; Cuuw 表 示 地 下 生物 碳 储量 (Mg); Coi ZEN 
土壤 矶 储量 (Mg) 5 Ci 表示 死亡 有 机 碳 储量 (Mg)。 
各 种 参数 通过 不 同文 献 搜 集 获取 ””, 具 体 结果 如 
表 2。 

本 研究 采用 InVEST 3.4 软件 (The Integrate Val- 
uation of Ecosystem Services and Tradeoffs Tool ) 进行 
碳 储量 的 估算 , 它 是 由 美国 斯 坦 福 大 学 ` 大 自然 保 
护 协 会 和 世界 自然 基金 会 联合 开发 的 开源 式 生 态 
系统 服务 功能 评估 模型 ,是 一 种 整合 自然 资本 作为 
决策 依据 的 先进 工具 ”” 。 

1.4 不 同 地 形 因子 对 碳 储量 的 影响 

坡度 、 坡 向 和 坡 位 通常 是 表述 地 形 因子 的 3 个 
主要 指标 。 麻 塔 流域 坡度 范围 为 0"~36.63" ,为 分 析 
不 同 坡度 范围 碳 储量 及 其 变化 ,将 坡度 分 为 0"~5"、 


5°~10°, 10°~15° ,15°~20° ,20°~25° ,.25°~30° 、>30° 等 
级 ;为 分 析 不 同 坡 向 碳 储 量 及 其 变化 ,将 坡 向 按照 
光照 由 弱 到 强 分 为 阴 坡 、 半 阴 坡 、 半 阳 坡 、 阳 坡 , 即 
流域 按照 自 正 北方 向 0°, 顺 时 针 旋 转 到 正 南方 为 
180°, 总 体 1 周 为 0"~360°, 其 中 0°~45° 与 315°~360? 
合 称 为 阴 坡 ,45°~90° 与 270°*~315° 合 称 为 半 阴 坡 ， 
90°~135° 与 225°~270° 合 称 为 半 阳 坡 ,135°~225° 合 称 
为 阳 坡 ;整个 流域 海拔 范围 为 1185~1379 m, 为 分 析 
不 同 坡 位 的 碳 储量 及 其 变化 ,依据 海拔 范围 等 间距 的 
将 坡 位 划分 为 3 个 梯度 ,其 中 海拔 1185~1250 m 称 为 
下 坡 位 海拔 1250~1315 m 称 为 中 坡 位 ,海拔 1315~ 
1379 m 称 为 上 坡 位 。 用 ArcMap 10.6 绘 制 麻 塔 流域 
碳 储量 空间 格局 图 ,用 及 软件 中 raster 包 分 析 不 同 地 
形 因子 对 碳 储量 的 影响 ,并 进行 相关 的 统计 作 图 。 


2 结果 与 分 析 


21 不 同时 期 流域 碳 储量 格局 及 其 变化 特征 
由 图 2 可 知 ,1999 年 生态 系统 总 碳 储 量 为 


R2 ， 麻 塔 流域 土地 利用 固 碳 参数 ( 碳 密度 ) 
Tab. 2 Carbon parameter for land use type in the Mata Watershed (carbon density) 


十 地 利用 类 型 编码 地 上 生物 碳 密度 地 下 生物 碳 密度 SBE I 死亡 有 机 碳 密度 
/(Mg+ hm’) /(Mg+ hm’) /(Mg:hm’) /(Mg+ hm’) 
森林 1 44.60 11.10 106.10 1.90 
灌木 2 7.59 4.00 80.50 0.80 
草地 3 1.34 10.07 85.40 0.45 
果园 4 14.10 5.60 76.00 0.03 
KH 5 3.00 0.68 76.83 0.07 
其 他 6 0.00 0.00 40.90 0.00 
(a) 1999 年 (c) 1999 一 2016 年 
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图 2 1999 年 和 2016 年 碳 储量 空间 格局 及 1999 一 2016 年 碳 储量 变化 
Fig. 2 Spatial pattern of carbon storage in Mata Watershed in the year of 1999 and 2016 and carbon 
storage change between 1999 and 2016 
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22120.25 Mg( 碳 密度 为 90.66 Mg-hm’), 2016 FÆ 
态 系统 总 碳 储量 为 23808.61 Mg ( 碳 密度 为 97.58 
Mg.hm)。 总 体 来 看 ,1999 一 2016 年 , 麻 塔 流域 总 碳 
储量 增加 1688.36 Mg( 碳 密度 增加 6.92 Mg hm’) ,总 
国 碳 功能 提升 约 7.63%。 流 域 生 态 系统 碳 储量 的 增 
加 主要 位 于 流域 的 南部 地 区 。 
2.2 不 同 土地 利用 类 型 的 碳 储量 及 其 变化 特征 
1999 一 2016 年 间 , 麻 塔 流 域 土地 利用 类 型 的 转 
换 导 致 各 土地 利用 的 碳 储量 变化 较为 剧烈 (图 3)。 
其 中 ,森林 .草地 .果园 碳 储量 有 增加 的 趋势 ,这 3 种 
土地 利用 类 型 的 碳 储量 总 体 提 升 约 7484.62 Mg ,其 
中 森林 占 总 提升 的 17.72% ,草地 占 总 提升 的 
42.57% ,果园 提升 占 总 体 提升 的 39.71% ; 而 灌木 AR 
田 和 其 他 土地 利用 从 1999 一 2016 年 的 碳 储量 减少 
T 5796.27 Mg, 其 中 灌木 减少 最 多 , 占 总 体 减 少 的 
76.01%。 这 表明 森林 .草地 和 果园 土地 利用 类 型 是 
1999—2016 年 麻 塔 流域 土地 利用 转换 后 流域 景观 


Dee EGET AY) EE TRA o 
12000 m 19994F 0 20164 
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#8 4000 
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森林 | 灌木 ”草地 农田 果园 其 他 
土地 利用 类 型 
图 3 1999 年 和 2016 年 麻 塔 流域 不 同 土地 利用 类 型 碳 储量 


Fig. 3 Carbon storage of different land use types in 1999 and 
2016 in Mata Watershed 


2.3 地 形 因子 对 碳 储量 及 其 变化 的 影响 

麻 塔 流域 7 个 坡度 等 级 梯度 碳 储量 及 其 变化 的 
结果 表明 (图 4a) ,1999 一 2016 年 该 流域 碳 储量 呈 先 
增加 后 减 小 的 趋势 , 旦 碳 储 量 主要 位 于 坡度 0°~25?° 
之 间 ( 占 整个 流域 总 碳 储量 的 90% 左 右 ) ,其 中 10°~ 
15" 坡 度 范 围 内 总 碳 储量 最 高 。 从 1999—2016 年 
间 , 碳 储量 增加 较 大 的 区 域 主要 位 于 10°~30° ,该 坡 
度 范围 内 占据 总 体 碳 储量 提升 比例 的 96.78% ,其 他 
坡度 在 过 去 提升 不 明显 或 略微 降低 (如 0°~5°), 这 
表明 坡度 10°~30° 是 土地 利用 结构 改造 和 固 碳 功能 
提升 的 主要 地 形 位 置 。 


研 F 
= 12000 (a) 坡度 199948: 
Pa 8000 201647: 
a 4000 1 
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图 4 麻 塔 流域 地 形 因子 对 碳 储量 的 影响 
Fig.4 Influence of topography factors on carbon storage in 
Mata Watershed 


化 的 结果 显示 (图 4b ) ,无 论 在 1999 年 还 是 2016 年 ， 
流域 生态 系统 碳 储量 主要 分 布 在 下 坡 位 和 中 坡 位 ， 
上 坡 位 占有 比例 较 小 。 在 不 同 坡 位 上 , 1999—2016 
年 碳 储量 提升 主要 在 中 坡 位 和 上 坡 位 (提升 了 
2354.07 Mg) , 而 下 坡 位 碳 储量 降低 (下 降 了 665.71 
人 结构 改造 
过 程 中 ,中 坡 位 和 上 坡 位 是 碳 储量 增长 的 主要 贡献 
地 形 部 位 。 

1999 一 2016 年 麻 塔 流域 总 碳 储量 在 不 同 坡 向 
的 分 布 格局 表明 (图 4c), 碳 储量 的 提升 主要 处 于 半 
阴 坡 ,而 其 他 坡 向 对 碳 储量 提升 较 小 或 基本 保持 稳 
定 ( 如 阳 坡 )。1999 一 2016 年 , 麻 塔 流域 半 阴 坡 提 升 
矶 储量 的 比例 占 总 体 提升 的 比例 的 70.71% ,表明 半 
阴 坡 是 土地 利用 转换 和 固 碳 功能 提升 的 主要 地 形 
位 置 。 


3 讨论 


3.1 1999—2016 年 麻 塔 流域 碳 储 量 的 变化 及 其 政 
策 解释 

本 研究 基于 RS 和 GIS 的 手段 ,利用 InVEST 模 
型 评估 了 1999—2016 年 麻 塔 流域 土地 利用 转换 后 
流域 整体 固 碳 能 力 及 其 变化 。 研 究 表明 :经 过 18 a 
的 土地 利用 改造 ,流域 整体 固 碳 功能 提升 约 
7.63%。 在 此 期 间 ,将 农田 、 灌 从 和 裸 地 土地 利用 类 
型 改造 为 森林 .草地 和 果园 土地 利用 类 型 是 流域 景 
观 碳 储量 提升 的 主要 直接 原因 ,尤其 是 流域 上 游 地 


麻 塔 流域 不 同 坡 位 的 碳 储量 分 布 格局 及 其 


区 (南部 ) 的 森林 面积 增加 和 流域 中 下 游 果 园 面积 
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的 增加 (图 1)。 

森林 .草地 和 果园 的 面积 增加 主要 发 生 在 坡度 
10°~30° .中 上 坡 位 的 半 阴 坡 , 从 而 直接 导致 了 流域 
整体 碳 储 量 的 增加 。 了 驱动 麻 塔 流域 土地 利用 变化 
的 社会 原因 主要 是 当地 的 退耕 还 林 (1999 年 开始 ) 
政策 和 近年 来 的 生态 农业 政策 。 在 植被 建设 过 程 
中 ,大量 的 农田 、 灌 从 和 裸 地 被 改造 为 林地 .草地 和 
果园 呈 。 本 研究 说 明了 退耕 还 林 工 程 过 程 中 将 农 
田 .灌木 和 裸 地 等 土地 转换 为 林地 有 利于 麻 塔 流域 
生态 系统 的 固 碳 , 更 重要 的 是 本 文 发 现 果园 也 是 麻 
塔 流 域 景观 碳 储量 提升 主要 植被 类 型 ,这 在 之 前 大 
量 的 关于 该 地 区 退耕 还 林 工 程 的 回 碳 效益 方面 研 
究 很 少 涉及 ”3 了。。 近 年 来 ,在 当前 新 农村 建设 与 生态 
农业 等 政策 的 驱动 下 ,将 延安 周边 地 区 的 水 土 保持 
林地 和 耕地 改造 为 果园 地 的 工作 逐渐 展开 。 本 研 
究 以 麻 塔 流域 为 案例 ,证 明了 退耕 还 果 政 策 与 退耕 
还 林政 策 相 似 , 也 有 利于 增强 该 地 区 流域 景观 的 固 
碳 效益 。 

虽然 退耕 还 林 ( 草 ) 和 退耕 还 果 政 策 能 够 驱动 
麻 塔 流域 景观 碳 储量 的 提升 ,但 提升 的 幅度 并 不 是 
很 高 ,这 主要 的 因为 以 麻 塔 流域 为 代表 的 黄土 丘陵 
区 为 干旱 、 半 干旱 区 ,生态 系统 潜在 固 碳 能 力 相 对 
于 其 他 湿润 地 区 较 低 。 例 如 麻 塔 流域 2016 年 的 平 
均 碳 密度 为 97.58 Mg.hm2( 根 据 流域 面积 和 碳 储量 
计算 得 到 ) , 低 于 广西 西江 流域 553(186.99 Mg-hm’) , 
贵州 晴 隆 号 (173.1 Mg hm?) 、 太 行 山 湛 河 流域 
(141.9 Mg hm’) .河北 黄 骅 “(112.34 Mg hm’) 等 
地 区 。 
3.2 不 确定 性 分 析 及 其 实践 意义 

本 研究 是 基于 InVEST 模 型 来 评估 1999—2016 
年 麻 塔 流域 碳 储量 及 其 变化 。 该 模型 内 在 的 假设 
认为 生态 系统 的 各 组 分 的 碳 密度 在 时 间 尺 度 上 是 
保持 不 变 的 小 。 随 着 植被 的 生长 ,生态 系统 各 组 分 
的 碳 密 度 可 能 会 持续 增加 。 同 时 ,InVEST 模 型 也 没 
有 考虑 人 对 农田 和 果园 的 管理 措施 可 能 会 提高 
态 系统 的 固 碳 潜力 。 因 此 ,评估 中 采用 固定 不 变 的 
矶 密度 值 可 能 会 低估 麻 塔 流域 碳 储量 的 潜力 。 

虽然 ,本 研究 并 不 能 准确 的 估算 麻 塔 流域 碳 储 
量 的 变化 ,但 我 们 认为 本 人 研究 能 够 反映 麻 塔 流域 过 
去 18 a 土地 利用 转换 后 碳 储量 的 变化 趋势 ,并 能 
帮助 我 们 在 空间 尺度 上 寻找 碳 储 量 增加 的 区 域 , 并 
获得 有 效 的 土地 利用 转换 策略 信息 。 例 如 流域 碳 


储量 在 坡度 、 坡 向 和 坡 位 梯度 上 的 分 析 表 明 ,坡度 
10°~30° , 半 阴 坡 和 中 上 坡 位 是 麻 塔 流域 碳 储 量 增长 
的 主要 贡献 地 形 部 位 。 这 意味 着 麻 塔 流域 在 这 些 
地 形 部 位 的 土地 利用 转换 策略 (转换 为 林地 或 果园 
地 ) 有 利于 生态 固 碳 效 益 的 提升 ,适合 于 广泛 的 
推广 。 

2016 年 的 麻 塔 流域 土地 利用 图 显示 这 些 地 形 
部 位 主要 种 植 了 果园 (以 苹果 为 主 ) ,经 济 果 园林 与 
水 土 保持 林 相 比 ,经 济 果园 林 固 碳 能 力 虽 弱 于 水 土 
保持 林地 (ES AS FE A AR Sh SAR 
田 、 草 地 和 灌木 地 被 转换 为 果园 是 麻 塔 流域 整体 碳 
储量 提升 的 主要 原因 。 

除了 考虑 其 生态 固 碳 效益 之 外 ,经 济 果园 林带 
来 的 经 济 效 益 不 容 忽视 ,只 有 保障 农民 经 济 收益 ， 
才能 使 水 土 保持 与 生态 修复 工程 更 加 长 远 和 持 
久 。 近 20 a 来 ,旱地 苹果 为 麻 塔 流域 居民 的 支柱 产 
业 , 近 年 来 扶贫 政策 的 提出 ,政府 对 麻 塔 流域 苹果 
产业 的 重视 ,形成 一 系列 完善 的 制度 和 政策 ,使 该 
流域 果园 面积 不 断 增 加 。 经 济 果林 面积 的 持续 增 
加 及 其 综合 优化 配置 ,将 是 未 来 黄土 丘陵 区 水 土 保 
持 及 经 济 发 展 的 趋势 ,也 必 将 成 为 黄土 高 原 未 来 发 
展 的 一 大 特色 。 增 强 该 地 区 经 济 果林 的 固 碳 效益 ， 
使 其 能 够 达到 水 土 保持 林 的 水 乎 ,使 生态 与 经 济 效 
益 最 大 化 ,未 来 还 有 更 多 的 研究 要 做 。 


4 结论 


(1) 麻 塔 流域 过 去 18 a 的 植被 结构 改造 中 ( 主 
要 是 苹果 林 建 设 ) ,区 域 植被 总 体 碳 储量 提升 了 
7.63% ,森林 .草地 和 果园 面积 的 急剧 增加 是 碳 储量 
增加 的 直接 原因 (其 中 果园 贡献 39.71% )。 

(2) 退耕 还 果 与 退耕 还 林政 策 对 麻 塔 流域 的 固 
碳 功能 的 影响 具有 相似 的 效益 ,有 利于 增强 该 地 区 
流域 景观 的 固 碳 能 

(3) 麻 塔 流域 碳 储 量 提升 主要 位 于 半 阴 坡 、 中 
上 坡 位 以 及 坡度 10°~30° ,说 明了 麻 塔 流域 在 该 地 
形 部 位 的 土地 利用 转换 模式 (转换 为 林地 或 果园 
地 ) 具 有 较 好 的 推广 价值 。 

(4) 经 济 果林 面积 的 持续 增加 及 其 综合 优化 配 
置 , 将 是 未 来 黄土 丘陵 区 水 土 保 持 及 经 济 发 展 的 趋 
势 , 增 强 该 地 区 经 济 果 林 的 固 碳 效 益 ,使 其 能 够 达 
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干旱 区 研究 
到 水 土 保持 林 的 水 平 将 能 够 使 该 地 区 的 生态 与 经 [11] Goldstein J H, Caldarone G, Duarte T K, et al. Integrating ecosys- 


济 效益 最 大 化 。 
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Abstract: The soil and water conservation forests in the areas surrounding Yan’an city are being replaced by 
economic forests such as orchard, which will definitely have an impact on the regional carbon storage. Therefore, 
it is necessary to determine a method to quantify the impact of orchard construction on the regional carbon 
storage capacity. Based on the INVEST model and ArcGIS platform, this study evaluated the change in carbon 
storage in the Mata Watershed in southern Yan’ an from 1999 to 2016, and explored the influence of slope, slope 
aspect, and slope position on the spatial distribution of change in carbon storage. The results showed that the 
carbon storage of forest, grassland, and orchard has increased, whereas that of shrub, cropland, and other lands 
has declined during the 18 years of vegetation and reconstruction in the Mata Watershed (mainly orchard 
construction). The regional total carbon storage has increased by 7.63%, of which 39.71% has been contributed 
by orchard. The increase in spatial position of carbon storage was mainly located in half-shaded-slopes, middle, 
and top of slopes; as well as slopes with range of 10°-30°. This study provides a methodological case study for 
quantifying the impact of economic forest construction on regional carbon storage function in the loess hilly 
region. 
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